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Abstract. - A lattice Boltzmann model is introduced which simulates oil-water-surfactant 
mixtures. The model is based on a Ginzburg-Landau free energy with two scalar order parame- 
ters. Diffusive and hydrodynamic transport is included. Results are presented showing how the 
surfactant diffuses to the oil-water interfaces thus lowering the surface tension and leading to 
spontaneous emulsification. The rate of emulsification depends on the viscosity of the ternary 
fluid. 



Introduction. - The addition of surfactant to a binary mixture of oil and water can produce 
many different complex structures on a mesoscopic length scale. The surfactant molecules move 
to the interface and lower the oil-water interfacial tension. This can result in, for example, 
lamellar, micellar, microemulsion or hexagonal arrangements of the oil and water domains ^. 

The equilibrium behaviour of such amphiphilic systems is well understood. However the 
dynamics of the self-assembly of the mesoscale phases and their rheology are less well investi- 
gated. This is a difficult problem because of the interplay between several relevant transport 
mechanisms, the diffusion of the constituent components and their hydrodynamic flow. To 
date models of amphiphilic rheology which treat hydrodynamic effects include time-dependent 
Ginzburg-Landau approaches^, ^, molecular dynamics and a lattice gas cellular automaton 
scheme based on microscopic interactions!^ 

The aim of this Letter is to introduce an alternative numerical scheme that can model the 
dynamics of amphiphilic systems in such a way that diffusive and hydrodynamic mechanisms 
are included. The numerical approach that we use is lattice Boltzmann simulations which have 
emerged as a useful tool to study the dynamics of complex fluids]^. We base our approach 
on that described by Orlandini et. al. ^] where the correct equilibrium of the fluid is 
imposed by choosing an appropriate free energy and including it in such a way that the fluid 
spontaneously reaches the equilibrium described by its minimum. 

Previous lattice Boltzmann models of amphiphilic systems have been based on a single order 
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parameter, that of the phase separating binary fluid |2|. The effect of the amphiphihc 
molecules has been mimiced by varying the surface tension in the input free energy. Although 
this approach proved successful it has the disadvantage of not including the surfactant dynam- 
ics explicitly. 

We first give a description of the method and then present results showing how the surface 
tension of the binary fluid interface is lowered by surfactant at a rate which depends on 
the surfactant diffusion constant. As the surface tension becomes negative, this leads to the 
break-up of the interface and to spontaneous emulsification JTsf to a lamellar phase. 

The Lattice- Boltzmann Scheme. - We consider a Ginzburg-Landau model defined by the 
free energy functional Q which depends on two scalar order parameters 0(r) and p(r) 

+Mp^ + p2<f>\V^p) + I33P^{V^)\ . (1) 

0(r) and p(r) can be identified, respectively, with the local density difference of oil and water 
and with the difference of local surfactant concentration from its average p. p enters the 
model via the parameter k. k is positive for small surfactant concentration. As k decreases 
and eventually becomes negative, p increases. 

The thermodynamic variables that we will need are the chemical potential difference be- 
tween oil and water A/Lt, the chemical potential A of the surfactant and the pressure tensor 
Pap- The chemical potentials follow from the free energy as [ p^ 

Am = ^ + 60^" kV20 + c(v20)2 + 2/3ip0 + 2/320(VV) + /33P(V20)+/33V2(p0), (2) 



A= — =ap~ AV^p + 7(V2p)2 + p^cj)^ + p^^^cj)^ + /Sg^lV^^). (3) 
dp 

The derivation of the pressure tensor is slightly more complicated. Considering a linear 
combination of all symmetric tensors having two or four gradient operators, we find that a 
suitable choice, which allows the pressure tensor to obey the equilibrium condition daPap ~ 
is 

Paf} - {pL + c[(V20)2 + a,^9,v20] -f 7[(V'p)' + ^d^pd^W^p] + f32 [d^^d^p + cj^^W^p] 

+P3[d„pq^d„<p + pq^v^4']} 

-l[d^pdpV^p+dppda,V^p] -f32[d^cj:'^dpp+d^pdpcj::^] - P^[dap4>dfj4>+da^dpp4>] (4) 

where 

PL = ^0^+ \b<^'- n^iV'^)~^iV^r+cHVr<t>- |(VV)^+ ^p^-XpiV'p) -^(Vpf 
+7P(V') V- 1 (V V)' + 2/3iP0' +/320' (V^p) -f /32p(v202) +;33P0(V20) (p0) (5) 

The lattice Boltzmann scheme is defined in terms of three distribution functions fi (r) , gi (r) 
and hi{r), each of which evolves during a time step At according to a single relaxation time 
Boltzmann equation |l^ 

f,(r + e,At,t + At)- Mr,t) = -l[/,(r, t) - /^(r, i)], (6) 
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g, {r + e,At,t + At)-g,{r,t) = -l[g^(r,t) - g°{r,t)], (7) 

h, {r + e,At,t + At) - h,{r,t) = - — [h,{r,t)-h°{r,t)] (8) 

where t, T0 and Tp are independent relaxation parameters and ei are the unit lattice vectors. 
The distribution functions are related to the physical variables by 

n = ^fi, nu^^fiBi, (/) = ^5i, p^^hi (9) 

i i i i 

where n is the total density and u is the mean fluid velocity. These quantities are locally 
conserved and, therefore, we require that the equilibrium distribution functions ff,g^,h}l also 
fulfil Eqs. (|^). The higher moments of /f, g^ and are defined by imposing the additional 
requirements 

f^eiaCip = Pq/3 + nUaUp , ^ g°eia = 4>Ua , ^ g°e.,aeif3 = Vcl,A^Ji5al3 + cj)UaUp , (10) 



/i°e„ = pua , ^ /i°eiQe,/3 = TpAiJa/s + (11) 

where and Fp are mobilities. The equilibrium distribution functions are defined as usual in 
terms of an expansion in the velocity u|^, We use a 9- velocity model on a square lattice 
to obtain the results presented here. 

These definitions lead to the continuum equations which follow from expanding Eqs. (^, 
(0) and d) toO(Ai2)0 

dtn -\- da{nua) = ^ dt{nua) + d(3{nUaUi3) — —dpPap + i^\7'^{nua) + da[X{n)dj{nUj)] , (12) 
dtp + do^ipu^) - Tp&pV^A - &pd^ (z^pPo^p) , (14) 



n 



where 

(2r-l),^^, , , 1,^^,1 dpo, ^ 1, ^ 1, 



V 



g (Ai), A(n) = (r - -)At{- - -^), 9^ = Ai(r^ - -), 9p = At(Tp - -). (15) 



Dynamical Behaviour. - Unless otherwise stated the simulations were run with a = — 1, 
5 = c = 1, K = 0.1, a = A = 7 = 1, /3i = 0, /32 -0.2, /Jg = 0.4, r = 100, = = 
(1 + l/\/3)/2 ||l and with units in which At = 1. 

Above the critical temperature (a > 0) the equilibrium configuration is a mixture of the 
three components, with constant n, (p — and p — 0. For small variations in (f), Eq. ( p^ ) 
can be linearised about (j) — resulting in a convection-diffusion equation with diffusion 
constant = aT^Q^. Similarly for small variations in p, Eq. ( p^ ) can be linearised resulting 
in a diffusion equation with Dp = aTpQp. Following we tested these predictions by 
measuring D^ and Dp as functions of T^Q^ and rp9p respectively by monitoring the decay 
of a sinusoidal perturbation. The measured values are shown in Fig. 1 for different values of 
a and a. Agreement with the predicted values is very good. 
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Fig. 1. - Diffusion constant D as a function of mobility TQ for temperatures above critical for (a) the 
oil-water density difference and (b) the surfactant density. The full lines are analytic results which 
follow from linearising Eqs. ( p^ ) and Data points were obtained by following the decay of a small 

sinusoidal perturbation. 
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Fig. 2. - Variation of (a) the surface tension of a flat interface, (b) the maximum surfactant density 
at the interface with time for different values of the surfactant mobility: • TpQp = 0, o FpOp — 0.1, 
★ TpQp = 2.5. The parameter k is fixed to the value —1.15 and F^Q^ — 0.1. 



Belovif the critical temperature (a < 0) the binary mixture phase separates into two distinct 
phases, symmetric about (f> — 0. Oil-water interfaces are now formed in the system. We next 
consider the diffusion of the surfactant towards these interfaces. A system set-up to contain 
an oil-water interface that was initially a hyperbolic tangent was initialised with p = and 
and then allowed to evolve. Figure 2a shows the surface tension plotted as a function of time 
for three different values of the surfactant mobility. For zero mobility the interface relaxes to 
its equilibrium shape thus causing a decrease in the surface tension. However, no surfactant 



A. LAMURA et al.: A LATTICE-BOLTZMANN MODEL OF TERNARY-FLUID MIXTURES 



5 




Fig. 3. Spontaneous cmulsification of a droplet as surfactant diffuses to the interface. Results 
arc given for the evolution with time (from top, t=46, 591, 1252, 2652, 5615) of (a) the oil-water 
density difference 4> for high viscosity; (b) the surfactant density p (grey-scahng from black white 
corresponds to minimum p => maximum p) for high viscosity; (c) 4> for low viscosity; (d) p for low 
viscosity, k = —1.15 and r,^G^ = ^pQp = 0.1. 

diffuses to the interface and the surface tension remains positive. For finite mobiUty however 
the surfactant can diffuse to the interface and the final surface tension is negative. The rate 
at which the final value is achieved clearly depends on the value of the diffusion constant. 
Similar considerations apply to the increase in the peak value of the surfactant density p at 
the interface as a function of time (Fig. 2b). 
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Once the surface tension has become negative the interface should become unstable [y_3[. 
Figure 3 shows the evolution with time of a drop of oil in water which is initially approximately 
circular. (Small randomness in the radius is necessary or the drop remains metastable.) The 
evolution of both the (j) and p fields are shown. Small perturbations on the surface of the drop 
grow to form lamellae of a width appropriate to that minimising the free energy. The budding 
of the drop initially reflects the symmetry of the underlying lattice. We believe this to be a 
consequence of the sharpness of the interface relative to the lattice spacing. 

Figure 3 compares the spontaneous emulsification at high (r = 100) and low (r = 0.585) 
viscosities. The change in viscosity makes a striking difference to the speed at which the 
lamellae form. 

Summary. - We have described a lattice Boltzmann model for ternary mixtures, such 
as oil-water-surfactant systems. The model includes hydrodynamic and diffusive modes and 
controls the fluid equilibrium via a chosen input free energy. The viscosity and diffusivities of 
the fluid phases can be changed thus allowing a study of their effects on the self-assembly and 
rheology of the amphiphilic fluid. In particular we presented results for the effect of surface 
diffusion on the rate of change of the interfacial tension and for the spontaneous emulsification 
which results as the tension becomes negative. 

The model is not trivial and much work remains to be done to explore its physical and 
numerical properties within a wide dynamic and static parameter space. A well-defined equi- 
librium and the ability to impose rather than measure transport coefficients are particularly 
useful in this task. Extensions to three dimensions and a closer comparison to experimental 
results are important endeavours. 

We thank Enzo Orlandini and Peter Coveney for helpful discussions. 
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